Load Packages and Data

Cleaning Data

Visualizations

buis_2022 %>% 
  mutate(species = case_when(
    sp == "HB" ~ "Eretmochelys imbricata",
    sp == "GR" ~ "Chelonia mydas")
  ) %>% 
  ggplot(aes(x = species, y = cclnuct, fill = species)) +
  geom_violin() +
  scale_fill_manual(values = c("chartreuse4", "hotpink4")) + 
  theme_minimal() +
  coord_flip () +
  labs(title = "Average Curved Carapace Length of Green and Hawksbill Sea Turtles",
       subtitle = "Reproductive Females Measured During 2022 Nesting Season on Buck Island, USVI",
       x = "Species", 
       y = "Curved Carapace Length (cm)",
       fill = "Species")
## Warning: Removed 10 rows containing non-finite values (`stat_ydensity()`).

buis_2022 %>% 
  mutate(species = case_when(
    sp == "HB" ~ "Eretmochelys imbricata",
    sp == "GR" ~ "Chelonia mydas")
  ) %>% 
  ggplot(aes(x = species, y = cclnuct, fill = species)) +
  geom_boxplot() +
  scale_fill_manual(values = c("chartreuse4", "hotpink4")) + 
  theme_minimal() +
  coord_flip () +
  labs(title = "Average Curved Carapace Length of Green and Hawksbill Sea Turtles",
       subtitle = "Reproductive Females Measured During 2022 Nesting Season on Buck Island, USVI",
       x = "Species", 
       y = "Curved Carapace Length (cm)",
       fill = "Species")
## Warning: Removed 10 rows containing non-finite values (`stat_boxplot()`).

GTBK_data %>%
  filter(especie_nombre_latino == "Chelonia mydas" |
         especie_nombre_latino == "Eretmochelys imbricata" |
         especie_nombre_latino == "Lepidochelys olivacea" |
         especie_nombre_latino == "Caretta caretta") %>%
  ggplot(aes(x = especie_nombre_latino, y = a_c_c, fill = especie_nombre_latino)) +
  geom_boxplot() +
  scale_fill_manual(values = c("goldenrod1","chartreuse4", "hotpink4","cornflowerblue")) +
  theme_minimal() +
  coord_flip() +
  facet_wrap(~ geografía) +
  labs(title = "Curved Carapace Length Measurements of Four Sea Turtle Species",
       subtitle = "Monitored by GTBK in Kino Bay, Mexico between 2010-2023",
       x = "Species",
       y = "Curved Carapace Length (cm)",
       fill = "Species")
## Warning: Removed 1 rows containing non-finite values (`stat_boxplot()`).

GTBK_data %>% 
  filter(numero_recaptura %in% c(2, 3, 4)) %>% 
  ggplot(aes(x = numero_recaptura)) +
  geom_bar()

GTBK_data %>% 
  filter(numero_recaptura %in% c(2, 3, 4)) %>% 
  ggplot(aes(x = numero_recaptura, fill = sexo_correcto)) +
  geom_bar()

GTBK_data %>% 
  filter(numero_recaptura %in% c(1, 2, 3, 4)) %>% 
  filter(especie_nombre_latino == "Chelonia mydas" | 
         especie_nombre_latino == "Eretmochelys imbricata" |
         especie_nombre_latino == "Caretta caretta" |
         especie_nombre_latino == "Lepidochelys olivacea") %>% 
  ggplot(aes(x = numero_recaptura, fill = especie_nombre_latino)) +
  geom_bar() +
  scale_fill_manual(values = c("firebrick2", "darkseagreen", "royalblue4", "yellow2")) +
  theme_minimal() +
  labs(title = "Number of Sea Turtle Recaptures by GTBK",
       subtitle = "Of Four Sea Turtle Species in Kino Bay, Mexico",
       x = "Number of Times Recaptured",
       y = "Number of Turtles",
       fill = "Species")

# recapture_tibble <- GTBK_data %>%
#   group_by(marca_nombre) %>%
#   select(nombre_tortuga, fecha, especie, area_de_monitoreo, marca_nombre)
# 
# 
# recapture_tibble %>%
#   group_by(marca_nombre) %>%
#   summarize(marca_nombre)
GTBK_data %>%
  group_by(especie_nombre_latino) %>%
  count (sexo_correcto)
## # A tibble: 9 × 3
## # Groups:   especie_nombre_latino [5]
##   especie_nombre_latino  sexo_correcto     n
##   <chr>                  <chr>         <int>
## 1 Caretta caretta        Indefinido        2
## 2 Chelonia mydas         Hembra          339
## 3 Chelonia mydas         Indefinido      483
## 4 Chelonia mydas         Macho            26
## 5 Eretmochelys imbricata Indefinido        7
## 6 Lepidochelys olivacea  Hembra           10
## 7 Lepidochelys olivacea  Indefinido        1
## 8 Lepidochelys olivacea  Macho             4
## 9 <NA>                   <NA>              8
# GTBK_2018_2023 %>% 
#   group_by(especie_nombre_latino) %>% 
#   count(Sexo)

# GTBK_2010_2018 %>% 
#   group_by(Especie) %>% 
#   count(Sexo) 

buis_2022 %>% 
  group_by(sp) %>% 
  count(sex)
## # A tibble: 2 × 3
## # Groups:   sp [2]
##   sp    sex        n
##   <chr> <chr>  <int>
## 1 GR    female    27
## 2 HB    female    29
GTBK_data %>%
  filter(numero_recaptura == 1) %>%
  group_by(especie_nombre_latino) %>%
  count (sexo_correcto)
## # A tibble: 9 × 3
## # Groups:   especie_nombre_latino [5]
##   especie_nombre_latino  sexo_correcto     n
##   <chr>                  <chr>         <int>
## 1 Caretta caretta        Indefinido        2
## 2 Chelonia mydas         Hembra          248
## 3 Chelonia mydas         Indefinido      380
## 4 Chelonia mydas         Macho            23
## 5 Eretmochelys imbricata Indefinido        7
## 6 Lepidochelys olivacea  Hembra           10
## 7 Lepidochelys olivacea  Indefinido        1
## 8 Lepidochelys olivacea  Macho             4
## 9 <NA>                   <NA>              1
GTBK_data %>%
filter(especie_nombre_latino != "NA") %>%
  filter(numero_recaptura == 1) %>%
ggplot(aes(x = sexo_correcto, 
                     fill = sexo_correcto)) +
  geom_bar() + 
  facet_wrap( ~ especie_nombre_latino, 
              # <- fct_relevel(especie_nombre_latino, 
              #                         "Caretta caretta",
              #                          "Eretmochelys imbricata", 
              #                          "Lepidochelys olivacea",
              #                          "Chelonia mydas"),
               scales = "free_y") +
   labs(title = "Species and Sex Distribution", 
       subtitle = "of turtles captured in Kino Bay", 
       x = "sex", 
       y = "number of individuals", 
       fill = "sex") +
 scale_fill_manual(values = c("#cd4071","#feca8d", "#721f81", "#000004"))

naniar::gg_miss_var(GTBK_data)

visdat::vis_miss(GTBK_data)
## Warning: Raster pixels are placed at uneven horizontal intervals and will be shifted
## ℹ Consider using `geom_tile()` instead.

visdat::vis_miss(GTBK_2018_2023)

Turtle Photos/Memes

A Green sea turtle hatchling headed toward the sea.

A Green sea turtle hatchling headed toward the sea.

A Hawksbill sea turtle laying her nest on Buck Island.

A Hawksbill sea turtle laying her nest on Buck Island.